Plastic flow in polycrystal states in a binary mixture 
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Using molecular dynamics simulation we examine dynamics in sheared polycrystal states in a 
binary mixture containing 10% larger particles in two dimensions. We find large stress fluctuations 
arising from sliding motions of the particles at the grain boundaries, which occur cooperatively to 
release the elastic energy stored. These dynamic processes are visualized with the aid of a sixfold 
angle aj(t) representing the local crystal orientation and a disorder variable Dj(t) representing a 
deviation from the hexagonal order for particle j. 

PACS numbers; 62.20.Fe,61.82.FU,61.43.-j 



Understanding the deformation mechanisms of poly- 
crystals under applied strain is of great scientific and 
technological importance in materials science [H, 0) S 0] ■ 
If the typical grain size exceeds a critical size dc ~ lOnm, 
interplay of grain boundaries and dislocation motions de- 
termines the mechanical properties. For very small grain 
sizes less than dc^ some simulations suggested that plas- 
tic deformations are caused by sliding motions at the 
grain boundaries. There remain a number of puzzles in 
polycrystal rheology not yet simulated microscopically, 
such as the Portevin-Le Chatelier effect in dilute alloys 
caused by intermittent yielding On the other hand, 
in physics, dynamics of sheared glassy materials has been 
studied extensively d, 0] , including microscopic particle 
systems [1, 0, , granular materials 11 , l3l , and foams 
[l3|, while not enough attention has been paid to poly- 
crystal rheology. 

Recently, using molecular dynamics simulation, we 
have examined the dynamics of polycrystal states real- 
ized on very small scales in a model binary mixture 14 1 . 
The parameters we have changed are the size ratio of the 
diameters of the two components (J2 j <y\ , the temperature 
T, and the composition c. Polycrystal states appear as 
intermediate states between crystal and glass, where the 
grain boundary motions are severely slowed down with 
size dispersity. The particles in the grain boundary re- 
gions have a relatively high mobility leading to dynamic 
heterogeneity on long time scales. Thus investigating the 
jamming dynamics over wide ranges of dij a\ and c helps 
us to understand the glass dynamics as highly frustrated 
limits. In this paper, we will present simulation results 
on polycrystal rheology, however, at fixed 02l(J\ and c. 

Our two-dimensional (2D) system is composed of the 
bulk region with volume V — L"^ and the top and bot- 
tom boundary regions with volume O.iy, as can be seen 
in Fig. 1. Shear flow was realized by the relative bound- 
ary motion of the boundaries, where the top and bottom 
velocities are ±7^/2, with 7 being an applied shear rate. 
In the bulk region, < ?/ < L, a mixture of large and 
small particles interact via a truncated Lenard-Jones po- 
tential of the form, Vafi(r) — A:e\{a j r)^"^ — (ca/j/r)^] — 
Caji ict,(3 — 1,2), characterized by the energy e and the 



soft-core diameter CTq^ — ((Tq -I- o'p)/'^ with (T2/(Ji = 1.4. 
For r > Tcut = 3.2(71, we set Vaf3{r) = and the con- 
stant Cap ensures the continuity of Vaf3{r) at the cut-off 
T = fcut- In this paper, the particle numbers are fixed 
at Ni = 900 and N2 = 100 or c = N2/{Ni + N2) = 0.1 
in the bulk. The volume V is chosen such that the vol- 
ume fraction of the soft-core regions is fixed at 1 or at 
4, = {Niaf + N2(tI)/V ^ \, so L ^ 33.1f7i. To each 
boundary (-0.1L <y<0 or L<y< I.IL), 100 small 
particles with radius ai are attached by the spring po- 
tential 10e|r — -Rjp. They also interact with the other 
particles in the boundary and bulk regions with the com- 
mon Lenard-Jones potential. Before our simulation, the 
attached positions Rj [j = 1, • • • , 100) in each bound- 
ary wall were determined in a liquid state realized at 
T = 2e/kB with the Lenard-Jones potential only. 

We inte grat ed the Newton equations using the leapfrog 
algorithm [l5| under the periodic boundary condition in 
the horizontal (y) direction, with the mass ratio being 
TO1/TO2 = (ci/(J2)^- The time step of integration is 
0.002r with 



= CTi y/mi/e. 



(1) 



We will measure the time t in units of r and the shear 
rate 7 in units of r^^. Without shear (7 = 0), (i) we first 
equilibrated the system in a liquid state at T = 2e/fcB in 
a time interval of 10'^ and then quenched it to the final 
temperature T — 0.2e/fcs. (ii) After a relaxation time of 
5 X 10^, there was no appreciable time evolution in various 
quantities obtained as an average over the particles (see 
Fig. 7 of Ref.Q). (ih) After these steps, we applied a 
constant shear to the system. 

Since a large fraction of the particles are enclosed by 
six particles in 2D dense particle systems, the local crys- 
talline order is represented by the sixfold orientation 
We define the orientation angle aj in the range [0,7r/3] 
for each particle j using the following complex number. 



*i = ^ exp[6i9jk] 



(2) 



/cGbondcd 



where the summation is over the particles "bonded" to 
the particle j. The two particles j d a and fc G /? are 
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bonded, if their distance \rj — r^l is shorter than Rap = 
1.25(Tct/3 0- The upper cut-oflF Rap is sUghtly longer than 
the first peak position of the pair-correlation function 
gapir)- The 9jk is the angle of the relative vector rj — rfc 
with respect to the x axis. Next we construct another 
nonnegative-definite variable representing the degree of 
disorder for each particle j by [IJI 



D, = 2 J2 [l-cos6K-a,)]- 



(3) 



fcGbondcd 



For a perfect crystal at low T this quantity arises from 
the thermal vibrations and is nearly zero, but for particles 
around defects it assumes large values in the range 5 — 20. 

In Fig. 1, we show the particle configuration with the 
large particles in blue, the angles aj, and the disorder 
variable Dj for sheared states with (a) 7 — 10^^, (b) 
10-3, and (c) 10-2. See Ref.Q for the color map of aj 
and Fig. 3 below for that of Dj. The polycrystal grains 
and the grain boundaries are unambiguously visualized 
using aj and Dj . By comparing the three kinds of snap- 
shots for each state, we notice that the large particles 
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FIG. 1; (Color on line) Particle configuration (left), sixfold 
orientation order (middle), and disorder variable (right) for 
three sheared states, where 7 = IQ-* (top), IQ-^ (middle), 
and 10"^ (bottom), in a binary mixture with c — 0.1 and 
Oil Ox = 1.4. The snapshots for 7 = IQ-^ and IQ-"^ are very 
similar, while for 7 = IQ-^ the displacement vectors of the 
particles in a time interval of 10 (right) show shear localization 
near the boundaries. 
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FIG. 2; (Color on line) Shear stress Cxyit) vs strain after 
application of shear in units of eaf^ with 7 = 10-^, 10""^, 
and IQ-^, exhibiting large temporal fluctuations and strong 
shear thinning. 



tend to form the grain boundaries. This tendency can 
be seen even without shear ^14] . but is intensified under 
shear. Remarkably, the grain structures are insensitive 
to 7 for not very large shear (7 < 10-'^ here), where the 
effect of the boundary walls does not extend into the bulk 
and the time-average of the horizontal velocity is linear 
with the gradient 7. To support this weakness of the 
boundary effect, almost the same grain structures were 
realized under the Lee-Edwards boundary condition [isj 
(not shown in this paper). However, for very large shear 
7 = 10-2 in (c), the velocity gradient becomes localized 
near the boundaries with larger crystalline regions in the 
middle being continuously rotated and deformed. 

In Fig. 2, the average shear stress <Jxy{t) is displayed in 
units of ecrj"2 as a function of the average strain 7t after 
application of shear at t = 0, which is the sum of the mi- 
croscopic shear-stress contributions over all the particles 
in the bulk divided by V 0, EBl ■ Each curve is a result 
of a single simulation run. It undergoes large tempo- 
ral fluctuations, arising from plastic deformations of the 
polycrystal structures. For not large shear, it varies from 
a minimum about 0.4 to a maximum about 1.0 on a time 
scale of order 0.1-7-^, and its time-average is insensitive 
to 7. The effective viscosity rjcs (time-average of (Jxy/i) 
thus behaves as 7^^. At the largest shear 7 — 10^^^ the 
time scale of the fluctuations is longer than 0.17-^ and 
(^xy{i) is larger than in the lower shear cases. 

In Fig. 3, we display Dj(t) at (a) t = 1110, (b) 1120, 
and (c) 1130 at 7 = 10-^, superimposing the displace- 
ment vectors Arj(t) = rj{t -f /Si) — rj{t) with At — 10. 
The color is given to each picture independently, ac- 
cording to its minimum and maximum of Dj. We also 
present cTxyit) in the corresponding time range. Between 
t = 1110 and 1120, the deformations are mostly "elastic" 
as in (a) and Oxyit) gradually increases on the average. 
However, between t = 1120 and 1130, the picture (b) in- 
dicates signiflcant " sliding" particle motions in the grain 
boundary regions, which are of order ui even for this 
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FIG. 3; (Color on line) Disorder variable Dj at (a) t = 1110, 
(b) 1120, and (c) 1130 at 7 = 10~^. The arrows represent the 
particle displacement Arj in the subsequent time interval of 
width 10, which are large in (b) for the particles in the grain 
boundary regions. The average stress (Jxy{t) is also shown in 
this time region. 



small At. These sliding motions are triggered collectively 
throughout the system (in our small system), leading to 
a catastrophic drop of axy{t). The particles writtten in 
orange, which are mostly larger ones, may be regarded 
to be in disordered configurations. Their number is of 
order 10 in (a) and (c), while it is about 50 in (b). In 
(c) noticeable particle displacements still continue. Large 
scale collective motions of the particles within the grains 
are also conspicuous, which was already noticed in our 
previous simulation without shear 

In summary, in polycrystal states with very small 
grains, we have found intermittent yielding on a time 
scale of 0.17"-"^ for not very large 7. It is caused by coop- 
erative sliding motions in the grain boundary regions in 
agreement with the atomistic simulations 
our small-scale simulation, however, we cannot determine 
the spatial scale of the cooperative sliding extending over 
grains, which should be relevant in real systems. With 
increasing c, the typical sizes of the crystalline regions 



become smaller and the stress fluctuations gradually de- 
crease. That is, weaker disorder results in larger stress 
fluctuations in plastic flow. The proportionality of the 
structural relaxation time to 7""'^ and the shear thinning 
behavior rjas oc 7"^ still hold for glasses, as observed ex- 
perimentally 17 1 and numerically [8|] . Similar plastic flow 
phenomena should be observable in colloidal mixtures on 
expanded scales, as a future experimental system. It is 
of great interest how our findings can be related to the 
other well-known examples of enhanced stress fluctua- 
tions, where the effect of the size dispercity should be 
further examined [1, [H, 12, l^ . 
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